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The interplay of charmonium production and suppression in In + In and Pb + Pb reactions at 
158 AGeV and in Au+Au reactions at ^/s = 200 GeV is investigated with the HSD transport 
approach within the 'hadronic comover model' and the 'QGP melting scenario'. The results for 
the J /*¥ suppression and the l P' to J/*¥ ratio are compared to the recent data of the NA50, NA60, 
and PHENIX Collaborations. We find that, at 158 AGeV, the comover absorption model performs 
better than the scenario of abrupt threshold melting. However, neither interaction with hadrons 
alone nor simple color screening satisfactory describes the data at y/s — 200 GeV. A deconfined 
phase is clearly reached at RHIC, but a theory having the relevant degrees of freedom in this 
regime (strongly interacting quarks/gluons) is needed to study its transport properties. 
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1. Introduction 

Measurements of charmonium production in heavy-ion collisions at different energies can 
provide clear signatures of the onset of deconfinement. Indeed, according to potential model pre- 
dictions and to the pioneering idea of Matsui and Satz [jl]], cc meson states might no longer be 
formed in a very hot fireball due to color screening This initially intuitive expectation has 

guided experimental studies for almost two decades. However, more recent lattice QCD calcula- 
tions have shown that the J /*¥ survives up to at least 1.5 T c (T c as 170 to 185 MeV) such that the 
lowest cc states may remain bound up to rather high energy density [||, ^, [7|, H]. On the other hand, 
the Xc an d appear to melt soon above T c . 

According to present knowledge, the charmonium production in heavy-ion collisions, i.e. cc 
pairs, occurs exclusively at the initial stage of the reaction in primary nucleon-nucleon collisions. 
At the very early stage color dipole states are expected to be formed (cf. Refs. |^,[T^]). These cc 
states are assumed to be absorbed in a 'pre-resonance state' before the final hidden charm mesons 
are formed. Such absorption - denoted by 'normal nuclear suppression' - is also present in p + A 
reactions and is determined by a dissociation cross section Ob ~ 4 to 7 mb. Those charmonia 
or 'pre-resonance' states that survive normal nuclear suppression during the short overlap phase 
of the Lorentz contracted nuclei furthermore suffer from (i) a possible dissociation in the decon- 
fined medium at sufficiently high energy density and (ii) the interactions with secondary hadrons 
(comovers) formed in a later stage of the nucleus-nucleus collision. 



In the QGP 'threshold scenario', e.g the geometrical Glauber model of Blaizot et al. [11] as 
well as the percolation model of Satz [Q], the QGP suppression '(i)' sets in rather abruptly as soon as 
the energy density exceeds a threshold value £ c , which is a free parameter. This version of the stan- 
dard approach is motivated by the idea that the charmonium dissociation rate is drastically larger in 
a quark-gluon-plasma (QGP) than in a hadronic medium [^]. On the other hand, the extra suppres- 



sion of charmonia in the high density phase of nucleus-nucleus collisions at SPS energies [O, 13 



14], [B|] has been attributed to inelastic comover scattering (cf. |£o|, [Hj, 0, |]| pi |(| |]1 H, H 
and Refs. therein) assuming that the corresponding //*P-hadron cross sections are in the order of a 



few mb [E4, 25J, |2g, 27Q . In these models 'comovers' are viewed not as asymptotic hadronic states 
in vacuum but rather as hadronic correlators (essentially of vector meson type) that might well sur- 
vive at energy densities above 1 GeV/fm 3 . Additionally, alternative absorption mechanisms might 
play a role, such as gluon scattering on color dipole states as suggested in Refs. [28, 29, 30, |T|] or 
charmonium dissociation in the strong color fields of overlapping strings [[32|]. 

We recall that apart from absorption or dissociation channels for charmonia also recombination 
channels such D + D —> X c + meson (X c = (7/*P,^ c , x P / )) play a role in the hadronic phase. These 
backward channels - relative to charmonium dissociation with comoving mesons - have been found 
to be practically negligible at the SPS energies [[33|], but extremely important at the top RHIC energy 



of = 200 GeV [p41- This is in accordance with independent studies in Refs. [26, 29, 35, 36] 



and earlier analysis within the HSD transport approach [ |37| , p8| ] . 

The explicit treatment of initial cc production by primary nucleon-nucleon collisions and the 
implementation of the comover model - involving a single matrix element Mq fixed by the data at 



SPS energies - as well as the QGP threshold scenario in HSD were explained in Ref. [33] (see Fig. 1 



of Ref. [p3] for the relevant cross sections). We recall that the 'threshold scenario' for charmonium 
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dissociation is implemented as follows: whenever the local energy density e(x) is above a threshold 
value Ej (where the index j stands for J /*¥ ,% c: m'), the charmonium is fully dissociated to c + c. 
The default threshold energy densities adopted are E\ = 16 GeV/fm 3 for £2 = 2 GeV/fm 3 for 
Xc, and £3=2 GeV/fm 3 for «F'. 

It is presently not clear, if also the D-mesons survive at temperatures T >T C , but strong corre- 
lations between a light quark (antiquark) and a charm antiquark (quark) are likely to persist 
One may also speculate that similar correlations survive also in the light quark sector above T c such 
that 'hadronic comovers' - most likely with different spectral functions - might show up also at 
energy densities above 1 GeV/fm 3 , which is taken as a characteristic scale for the critical energy 
density. Therefore, we study both possibilities: with and without comover absorption (and D + D 
recombination) at energy densities above the cut-energy density parameter e cut = 1 GeV/fm 3 . 

Since we aimed to answer, whether the charmonium dissociation mechanism is identical at 



SPS and top RHIC energies, we adopted in [ |34| ] the same cross sections for the color-dipole disso- 
ciation with nucleons as well the dissociation cross sections with comovers as in Ref. [|3|] for SPS. 
Consequently no free parameters entered our studies at the RHIC energy. We note that the hadronic 
comover reactions for the recreation of charmonia 7/*P, % c , *P' by D + D reactions are incorporated 
in all simulations. This is a 'default' in the comover absorption and recreation scenario and 'nec- 
essary' in the QGP 'threshold scenario' because (in view of Fig. ||, l.h.s.) practically all charmonia 
are dissolved due to the very high initial energy densities. Therefore, any model without recreation 
of charmonia is clearly ruled out by the PHENIX data. 

2. Comparison to data 

We directly step on with results for the charmonium suppression at SPS energies in comparison 
with the experimental data from the NA50 and NA60 Collaborations. These Collaborations present 
their results on 7/*P suppression as the ratio of the dimuon decay of 7/*P relative to the Drell- 
Yan background from 2.9 - 4.5 GeV invariant mass as a function of the transverse energy Ej, or 
alternative, as a function of the number of participants V part , i.e. 

B^CJ{J/*¥)/g(DY)\2,9-4.5, (2.1) 

where B^^ is the branching ratio for J /*¥ — > pL + pL~ . In order to compare our calculated results to 
experimental data, we need an extra input, i.e. the normalization factor Bu.u.Gnn{J pV) / 'Gnn(DY), 
which defines the 7/*P over Drell-Yan ratio for elementary nucleon-nucleon collisions. We choose 



B^O NN (J 1^) / o NN (DY) = 36 in line with the NA60 compilation Q15Q. 



Furthermore, the \P' suppression is presented experimentally by the ratio 

(2.2) 



B^'^nn)a^")/a(DY] 



B^(J/V^w)o(J/y)/o(DY)- 

In our calculations we adopt this ratio to be 0.0165 for nucleon-nucleon collisions, which is again 
based on the average over pp,pd,pA reactions [|42|]. 

We first show in Fig. [T] the calculated ratio ( |Q| ) as a function of N pan for Pb+Pb and In+In 
collisions at 158 A-GeV (upper plots) in the nuclear suppression scenario, i.e. without comover 
dissociation or 'QGP threshold suppression'. The dashed (blue) lines stand for the HSD result 
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Figure 1: The ratio Bjxu.a{J /a(DY) as a function of the number of participants in In+In (l.h.s.) and 
Pb+Pb reactions (r.h.s.) at 158 A-GeV. The full symbols denote the data from the NA50 and NA60 Collabo- 
rations (from Refs. [|40|, |i"5, 41 1), while the dashed (blue) lines represent the HSD calculations including only 
dissociation channels with nucleons. The lower parts of the figure show the HSD results in the same limit 
for the M*' to J '/*}? ratio as a function of N part (for In+In) or the transverse energy Et (for Pb+Pb). The solid 
(red) lines show the HSD results for the comover absorption model with a matrix element squared |Mo| 2 
= 0.18 fm 2 /GeV 2 . The (light blue) bands in the upper parts of the figure give the estimate for the normal 
nuclear y/^P absorption as calculated by the NA60 Collaboration. The vertical lines on the graphs reflect the 



theoretical uncertainty due to limited statistics of the calculations. The figure is taken from [33 1 



while the (light blue) bands give the estimate for the normal nuclear absorption as calculated 
by the NA60 Collaboration. The normal nuclear suppression from HSD is seen to be slightly lower 
than the (model dependent) estimate from NA60, however, agrees quite well with their model 
calculations for more central reactions. The various experimental data points have been taken 

As a next step we add the comover dissociation channels within the 

0.18 fm 2 /GeV 2 . Note that in this 



from Refs. |pL HOL Hit]. 



model described in [33] for a matrix element squared |Mo| 2 
case the charmonium reformation channels are incorporated, too, but could be discarded since the 
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Figure 2: Same as Fig. |T] but for the 'QGP threshold scenario' with Ej/y = 16 GeV/fm 3 , e Zc = 2 GeV/fm 3 
= &p/ while discarding comover absorption. The figure is taken from p3|]. 



charmonium regeneration is negligible at SPS energies (cf. Ref. [0). The extra suppression of 
charmonia by comovers is seen in Fig. [j] (solid red lines) to match the J/^V suppression in In+In and 
Pb+Pb as well as the to //¥ ratio (for Pb+Pb) rather well. The more recent data (1998-2000) 
for the *+" to J /*¥ ratio agree with the HSD prediction within error bars. This had been a problem 
in the past when comparing to the 1997 data (dark green stars). The *+" to J /*¥ ratio for In+In 
versus centrality is not yet available from the experimental side but the theoretical predictions are 
provided in Fig. [I] and might be approved/falsified in near future. 

The results for the 'threshold scenario' are displayed in Fig. § in comparison to the same 
data for the thresholds Ejw = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 = e^. In this scenario the J/m 
suppression is well described for In+In but the suppression is slightly too weak for very central 
Pb+Pb reactions. This result emerges since practically all % c and \F' dissolve for N part > 100 in 
both systems whereas the J ' p¥ itself survives at the energy densities reached in the collision. Since 
the nucleon dissociation is a flat function of N part for central reactions, the total absorption strength 
is flat, too. The deviations seen in Fig. |2| might indicate a partial melting of the J /*P for N part > 
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Figure 3: The ratio B fJfi O'(.// v F) /a{DY) as a function of the number of participants N par t in In+In (red line 
with open squares) and Pb+Pb reactions (blue line with open circles) at 158 A-GeV relative to the normal 
nuclear absorption given by the straight black line. The full dots and squares denote the respective data from 
the NA50 and NA60 Collaborations. The model calculations reflect the comover absorption model (right 
part) and the 'QGP threshold scenario' (left part) with Ejpp = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , %/ = 6.55 
GeV/fm 3 while discarding comover absorption. Figure is taken from [^3||. 



250, which is not in line with most lattice QCD calculations claiming at least Ej/xy > 5 GeV/fm 3 . 
In fact, a lower threshold of 5 GeV/fm 3 (instead of 16 GeV/fm 3 ) for the J /*¥ has practically no 
effect on the results shown in Fig. ||. Furthermore, a threshold energy density of 2 GeV/fm 3 for 
the *F' leads to a dramatic reduction of the to ratio which is in severe conflict with the data 
(lower part of Fig. 0). Also note that there is no step in the suppression of J versus centrality. 



As pointed out before by Gorenstein et al. in Ref. [43], this is due to energy density fluctuations in 
reactions with fixed N part (or Et). 

Additionally, one can plot the results in an intuitive though model-dependent way, as a ratio of 
the measured J f*¥ yield divided by the normal nuclear absorption result calculated in the Glauber 
model. Since the NA60 Collaboration prefers to represent their data in this form, we additionally 
show in Fig. ^| our calculations for In+In (red lines with open squares) and Pb+Pb (blue lines 
with open circles) as a function of the number of participants N part relative to the normal nuclear 
absorption given by the straight black line 1 . The full dots and squares denote the respective data 
from the NA50 and NA60 Collaborations. The model calculations reflect the comover absorption 
model (right part) and the 'QGP threshold scenario' (left part) with Ej/y = 16 GeV/fm 3 , e Xc = 
2 GeV/fm 3 , ey = 6.55 GeV/fm 3 . Since only the representation is different the message stays the 
same: The comover absorption model follows slightly better the fall of the J /*¥ survival probability 
with increasing centrality whereas the 'threshold scenario' leads to an approximate plateau in both 
reactions for high centrality. 

Let us now move to a much higher energy scale by calculating charmonium dynamics at the 



1 Note that recently the NA60 collaboration has refitted the parameters of their Galuber model, therefore newer data 
releases might appear to be up- or down-scaled compared to the data plotted here jl5|], if shown in this particular 
representation (measured to expected ratio). This scaling falls within the systematic uncertainty of the ratio and does not 
change results and conclusions of our study. 
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Figure 4: The J/*¥ nuclear modification factor Raa for Am +Au collisions at y/s = 200 GeV as a function 
of the number of participants N par t in comparison to the data from [10] for midrapidity (full circles) and 
forward rapidity (full triangles). HSD results for the QGP 'threshold melting' scenarios are displayed in 
terms of the lower (green solid) lines for midrapidity J/*¥'s (\y\ < 0.35) and in terms of the upper (orange 
dashed) lines for forward rapidity (1.2 < y < 2.2) within different recombination scenarios (see text). The 
error bars on the theoretical results indicate the statistical uncertainty due to the finite number of events in 
the HSD calculations. Predictions for the ratio (I*') ay / B ^{J / l P)a> />j> as a function of the number of 
participants N par t are shown in the lower set of plots. The figure is taken from [p4Q. 



top RHIC energy of y/s = 200 GeV. In the initial stages of Au + Au collisions at this y/s, energy 
densities above 30 GeV/fm are reached [[34]]. Therefore, in the threshold melting scenario, all 
initially created J pif, and % c mesons melt. However, the PHENIX collaboration has found that 
at least 20% of J do survive at RHIC [45]. Thus, the importance of charmonium recreation 
is shown again. We account for J /*¥ recreation via the DD annihilation processes as explained 
in detail in [ |33| , [PJ. Note that in our approach, the cross sections of charmonium recreation in 
D + D — > / p¥ + meson processes is fixed by detailed balance from the comover absorption cross 
section J + meson — > D + D. But even after both these processes are added to the threshold 
melting mechanism, the centrality dependence of the Raa(J cannot be reproduced, especially 
in the peripheral collisions (see Fig. |j). This holds for both possibilities: with (r.h.s. of Fig. |]) and 
without (center of Fig. |]) the energy density cut £ cut , below which D-mesons and comovers exist 
and can participate in D + D <-> + meson reactions. 
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Figure 5: Same as Fig. 4 for the 'comover absorption scenario' including the charmonium reformation 
channels without cut in the energy density (l.h.s.) and with a cut in the energy density e cut = 1 GeV/fm 3 (see 
text for details). The figure is taken from [p4|. 



We recall that the nuclear modification factor Raa is given by 

KAA N colr dN(J/V) pp /dy' ^ 

where dN(J /*¥)AA/dy denotes the final yield of J /m in AA collisions, dN(J /^pp/dy is the yield 
in elementary pp reactions, N co u is the number of binary collisions. 

Comover absorption scenarios give generally a correct dependence of the yield on the central- 
ity. If an existence of D-mesons at energy densities above 1 GeV/fm 3 is assumed, the amplitude of 
suppression of J /m at mid-rapidity is also well reproduced (see the line for 'comover without £ CHf ' 
scenario in Fig.|5j, l.h.s.). Note that this line correspond to the prediction made in the HSD approach 



in [|38|]. On the other hand, the rapidity dependence of the comover result is wrong, both with and 
without e cut . If hadronic correlators exist only at e < e cut , comover absorption is insufficient to 
reproduce the J /*¥ suppression even at mid-rapidity (see Fig. ||, r.h.s.). The difference between the 
theoretical curves marked 'comover + £ cuf ' and the data shows the maximum supression that can 
be attributed to a deconfined medium. 
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3. Summary 

We have investigated the formation and suppression dynamics of J/*¥, % c and *F' mesons 
within the HSD transport approach for In + In and Pb + Pb reactions at 158 AGeV and for Au + Au 
reactions at y/s = 200 GeV. Two currently discussed models, i.e. the 'hadronic comover absorption 
and reformation' model as well as the 'QGP threshold melting scenario' have been compared to the 
available experimental data. We adopted the same parameters for cross sections (matrix elements) 
or threshold energies at both bombarding energies. 

We find that both scenarios are compatible with experimental observation of J/W suppression 
at SPS energies, while the *F' to J /*¥ ratio data appear to be in conflict with the 'threshold melting' 



scenario Q33|]. On the other hand, both 'comover absorption' and 'threshold melting' fail severely 
at RHIC energies [p4[|. The failure of the 'hadronic comover absorption' model goes in line with 
its underestimation of the collective flow V2 of leptons from open charm decay as investigated in 



Ref. [|46p. This suggests that 1) a deconfmed phase is clearly reached at RHIC, 2) the dynamics of 
c,c quarks at this energy are dominated by partonic interactions in the strong QGP (sQGP) which 
cannot be modeled by 'hadronic' interactions or described appropriately by color screening alone. 
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